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We study the stability under linear perturbations of a class of static solutions of 
Einstein-Gauss-Bonnet gravity in D = n+2 dimensions with spatial slices of the form 
E™ x M + , E™ an n— manifold of constant curvature k. Linear perturbations for this 



class of space-times can be generally classified into tensor, vector and scalar types. 



In a previous paper, tensor perturbations were analyzed. In this paper we study 
vector and scalar perturbations. We show that vector perturbations can be analyzed 
in general using an S-deformation approach and do not introduce instabilities. On 
the other hand, we show by analyzing an explicit example that, contrary to what 



happens in Einstein gravity, scalar perturbations may lead to instabilities in black 

holes with spherical horizons when the Gauss-Bonnet string corrections are taken 

CN ■ into account. 
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I. INTRODUCTION 

The analysis of the properties and behavior of gravity in higher dimensions has become 
a major area of research in recent years, motivated in particular by developments in string 
theory. Among others, the Einstein- Gauss-Bonnet (EGB) gravity theory has been singled 
out as relevant to the low energy string limit [1]. The EGB theory is a special case of Love- 
lock's theory of gravitation [2], whose lagrangian is a linear combination of Euler densities 
continued from lower dimensions. Lovelock's theory gives equations involving up to second 
order derivatives of the metric, and has the same degrees of freedom as ordinary Einstein 
theory [2]. A number of solutions to the EGB equations, many of them relevant to the 
development of the AdS — CFT correspondence [3], are known, among them a variety of 
black holes in asymptotically Euclidean or (A)dS spacetimes [4-9]. These were found mostly 
because they are highly symmetric. Analyzing their linear stability, however, confronts us 
with the complexity of the EGB equations, since the perturbative terms break the simplify- 
ing symmetries of the background metric. To be more specific, we consider spacetimes that 
admit locally a metric of the form 

ds 2 = -f{r)dt 2 + g{r)dr 2 + r 2 g lj dx i dx j , (1) 

where gijdx % dx^ is the line element of an n— dimensional manifold of constant curvature 
k — 1,0 or —1. Linear perturbations around (1) can be conveniently classified, following 
the scheme proposed in [10], into tensor, vector, and scalar perturbations. The k = 1 case 
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X™ = S n gives, for appropriate / and g, cosmo logical solutions, as well as higher dimen- 
sional Schwarzchild black holes. The stability of these solutions under tensor perturbations 
was studied in [11]. A complete classification of solutions to the EGB equation with line 
element (1) and k — 1,0 or — 1, together with a case by case stability study under tensor 
perturbations appears in [12]. In this paper we extend the analysis in [12] to vector and 
scalar perturbations. The methods we employ here are rather different from the analytic 
computations carried out in [11] and [12] for arbitrary dimension. The complexity of the 
computations involved in the study of vector and scalar perturbations forced us to develop 
alternative approaches. Explicit expressions were worked out for space-time dimensions 
D = n + 2 < 11, and their n— dependence was then interpolated. As a result, we obtained 
expressions which are low degree (< 4) polynomials in n, correct at least up to the highest 
physically interesting space-time dimension (D = 11, n = 9). Nevertheless, we conjecture 
that they are correct for arbitrary n. This is because only the lower values of n (up to n = 6) 
are required to obtain the n— dependence, and higher values appear as "predictions". More- 
over, it is clear that a purely analytic derivation of the equations for general n should lead 
to expressions of the form presented here, and this gives further support to our conjecture. 
The paper is organized as follows. In section II we give the general solution of the EGB 
equations with metric (1), then we introduce vector and scalar perturbations following [10]. 
In section III we study the stability of (1) under vector perturbations in EGB gravity. As 
was done for tensor perturbations in [12], we show that the problem reduces to obtaining 
a lower bound for the spectrum of a Schrodinger operator which, in spite of its complexity, 
turns out to be an S— deformation [10] of a much simpler operator. Section IV is devoted 
to scalar perturbations. A step by step guide for constructing the potential of the associ- 
ated Schrodinger operator is given, since its general expression for arbitrary dimensions is 
extremely long and has a complicated dependence on the parameters. This makes a general 
analysis of the stability problem under scalar type perturbations practically impossible. On 
the other hand, since the procedure given in this paper allows for the explicit construction 
of the potential for any particular choice of dimension n, and other parameters of the theory, 
our results can easily be used in the analysis of particular classes of situations. As an ex- 
ample and application, we analyze Schwarzschild-like five dimensional black holes and find 
a scalar instability of low mass, which is absent in Einstein gravity. Section V contains the 
conclusions. 

II. VECTOR AND SCALAR PERTURBATIONS OF A CLASS OF STATIC 

SPACETIMES 

The Einstein- Gauss-Bonnet (EGB) equations in the vacuum case are given by 
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Here we consider static spacetimes satisfying (2) with a metric of the form (1). These 
are foliated by spacelike hypersurfaces, orthogonal to the time-like Killing vector d/dt, that 
contain a submanifold of dimension n = D — 2 (D the spacetime dimension) of constant 
curvature k — 1,0 or —1, and line element g^dx % dx d . The non-zero components of the 
Riemann and Ricci tensors, and the Ricci scalar for a metric of the form (1) are given in 
[12]. Inserting these in (2) we find that (1) solves the EGB equation (2) if [4, 5, 7, 8, 11, 12] 

_L = f{r) = K _ r ^{r) (4) 
g{r) 

and ip(r) is a solution of 

4 y 2a (n + l)a ar n + l v ; 



This implies that, 



j, / r 2 / I 4a (n- I n- 2 

f ( r ) = K + _ 1 + e J H i — - — ^— - 

17 V ; a(n- l)(ra-2) I V n(n + l) 



(6) 



V> < and > PM,) (if «=-l) (10) 



where e = ±1 picks a root of the quadratic polynomial P{jp). The Ricci scalar for this 
solution is [12] 

R = (n + 2)(n + l)^(r) + 2r(n + 2)^> + r 2 ^!^, ( 7 ) 

ar dr z 

so that there is a singularity whenever ^ — > oo or ^ — > the stationary point of P, since 
dip /dr — > oo at this point. We also note [12] that if /x/a > the condition f — k — r 2 ip > 
reduces to 

V> < or < V, < PM) (if k=1) (8) 

a 

^ < ( if k = 0) (9) 

whereas for /i/a < 0, / = k — r 2 ip > is equivalent to 

ip < or < ip, P(ip) < -M*^ (if k = 1) (11) 

a 

if> < (if re = 0) (12) 

i) < and P(^)>-H^ (if re = -l) (13) 

a 

In [12], the space-times (l)-(4)-(5)-(6) were classified by studying the intersections of the 

ro+1 

curves P{ip) and ■ 

To study linear perturbations we follow the treatment given in [10], where it is shown 
that an arbitrary perturbation of the metric is a linear combination of perturbations of the 
tensor, vector and scalar types. Tensor perturbations around the EGB vacuum solution 
(l)-(4)-(5)-(6) were studied in [11, 12]. We consider now vector and scalar perturbations. 
We use a, b, c, d, ... as generic indices, greek indices /i, v refer to r, t, whereas i, j, k, I, m, ... are 
assumed to take values on A bar denotes tensors and operators on The perturbations 
are of the form 

9ab -> 9ab + K b , (14) 
and indices of h ab are raised using the background metric, therefore 5g ab = —h ab . 
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A. Vector perturbations 

A general vector perturbation is given by, 

V = h lti = rf,y i , h ij = 2r 2 H T V ij (15) 

where / M , and H T are functions of (r, t), and Vi, V^- are defined on S™. is a divergence-free 
vector harmonic field, satisfying, 

(A + k 2 y)Vi = 

ViF* = (16) 

where A := VjV 7 and V« are, respectively, the Laplacian and covariant derivative for the 
metric of (1). The symmetric tensor 

V ij = -^-(ViV j + V j Vd (17) 

is a harmonic tensor on E™, with the properties, 

[A + 4 - (n + 1)k] ^ = (18) 

n = , = (19) 

As shown in [10], for fcy 7^ (n — the combinations, 

Fp = U+ ^rd,H T (20) 
are a basis for gauge invariant variables. 

B. Scalar perturbations 

Scalar perturbations are of the form 

V = F^S, V = r:F n S ii h ij = 2r 2 (i/^ cjijS + H T Sij). (21) 
In (21), S is a scalar harmonic 

(A + k 2 s )S = 0, (22) 
and one constructs scalar-type harmonic vectors and tensors 

Si = ~%S, S :j ^V;Y,S ■ -g tj S, (23) 

k s k s n 

which satisfy 

[A + kg - (n - 1)k] Si = 0, ViS* = ST.* = 0, (24) 

and 

V,-# = ( n ~ ~ nK) ^ [A + k 2 s - 2uk] S l3 = 0. (25) 
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Harmonic symmetric tensors of arbitrary rank on n— spheres are constructed in [13]. In 
general, we lack explicit expressions for harmonic tensors on arbitrary constant curvature 
manifolds, although bounds on the Laplacian spectrum can be derived. 

It can be shown [10] that, in analogy, and, extending to higher dimensions the well known 
Regge- Wheeler gauge for "even" perturbations in four space-time dimensions, one can, by 
an appropriate coordinate transformation, choose a gauge where T a = 0, and Ht = 0. This 
makes T a b and Hl a basis for gauge invariant quantities. We shall make this choice of gauge 
in what follows, but, for convenience, we rewrite T a b in the form, 



r 



dR 



rt 



rt 



dt 



Fit — fF tt 



(26) 



III. VECTOR PERTURBATIONS IN EGB GRAVITY 



In this section we study vector perturbations of (1) as defined in equations (15)-(20). As 
is clear from the derivations in [10], one can always choose a gauge where Ht = 0. This 
generalizes the Regge- Wheeler gauge for odd perturbations to higher dimensions, but more 
important, it simplifies considerably the analysis. With this simplifying gauge choice, after 
a long computation we find that the non-trivial linearized Einstein-Gauss-Bonnet equations 
5Gij = are equivalent to, 



dr 



F r j(r^+a(n-2)±(r^(K-f))) 
~ [j (r (n - 2) + a (n - 2) ± (r^ (« - /) 



(27) 



while 5Gi. r = can be written as 



d 2 F t d 2 F, 



r — 



dt dr dt 2 



1 " ; 'r 2 • n(n !)(// 2)(/,- /)) 



r - 2 



dt 



+F r f 



[n - 1) (n - 2) r 2 (« - /) - 2 (n - 1) r 



df 4 d 2 f 



dr 



dr 2 



+ (n- l)(n-2) 



1 



+§("-3) (n-4) (/ - K ) 2 



a 



+ 



n(„-2) ; (// :',)(/ /,) ■ r'j-\ r 2 



(k 2 v -(n-l)K) -2r 4 A| 



(28) 



and 5G it = leads to, 
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d 2 f 



dr 2 



+(n - l)(n - 2) 
+ (n - 3) 



+ (n-l) (^2r^+n(/-/t) + 2 

rf/ (/-«)) +2(n-2)r|/ 
a - 2 r 4 A 
(4 - (n - 1)«) 



dr V dr 



|(.-/) 2 -2-2r| + 2/ 



a(n-2) (r|J - ( n - 3) (« - /)) 

+/ r(4_n) { r2 |: r( ™" 4) + a (n - 1) (n - 2) (« - /)) 

(r 2 + a(n- 1) (n - 2) (k - /)) 



9 /9F, 



<9r \ <9t 



r r (n-2) ^2 



(29) 



As explained in the Introduction, the n dependence in these formulas was obtained by inter- 
polating the results for different n values, and checked in the physically relevant range n < 9 
and isolated higher n values. No more than four points were required for the interpolation 
in every case. We conjecture that these formulas are valid for every n. As required for 
consistency, (27), (28), and (29) are not independent. Setting 



Ft=p[r)f(r) — — — , F r 



dr 



P(r) dx(t,r) 
f(r) dt ' 



with 



p(r) 



d 



>" 2 ) - a {n -2)j- [r (n ' 3 \f(r) - «)] 



solves trivially (27) and makes both (28) and (29) equivalent to 



<9 2 $(t,r) r2 d 2 <S>(t,r) „ . s d$(t,r) , , ft 

- / Q _, + Qv(r) V ' + gy(r)$(t, r) = 0. 



dt 2 



dr 2 



Here, 



and 



where 



H 



$(t,r) 



dr 

dx(t,r) 
dt 



Qv(r) = 



f(k 2 v -(n-l)K)H 



a 2 n (n — 2) 2 (n-3) (n- l)^ 2 + 2 an (n — 2) (n — 3) ip + 4 (n-2)A« + 2n 

2n((n-l)(n-2)^a + l) 2 



(30) 
(31) 

(32) 

(33) 
(34) 

(35) 
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Also 

Qv(r) 



f 



df 



a (n — 2)r ( (n — 3)(/c — /) — r— ) + r 3 



x < a(n - 2) 



2f^ 2 / 2 



<ir 2 



— r 



df_ 

dr 



_/( n _2)( n -3)(«-/) 



+ (^)((ra-l)/ + ra-3)r 



, 3 rf / r 2 

+ r jnr 

dr 



If we introduce a Regge- Wheeler "tortoise" coordinate r*, such that, 

dr* 1 



dr /(?")' 

an "integrating" factor K(r(r*)), and we also separate variables 

1 



(36) 



(37) 



$(t,r*) = 



K v (r(r*)) 



0(r(r*))e^, 



(38) 



we find that the choice 

K v (r) 



r n + a(n - 2)r 2 -^-(r^ n - 3 \K - /)) 
or 



-1/2 



(39) 



reduces (32) to a (stationary) Schrodinger equation, 



H(f> = 



d 2 



+ W0 = -uJ 2 (j) = E<f> 



dr 



*2 



with "energy" eigenvalue — u 2 and "potential" V(r(r*)) given by 



V v (r) = q v {r) + 



Kvtr) 



d 2 K v 
dr 2 



dK v \ df 
dr } dr 



(40) 



(41) 



The stability problem therefore reduces to analyzing the spectrum of the Hamiltonian op- 
erator (40), the presence of a negative eigenvalue (oo real) signaling an instability. H acts 
on square integrable functions on the / > region /, and, in spite of the complexity of 
the potential Vy, information on its spectrum can be obtained by using the S— deformation 
approach, as done in [10, 12]. As in [12], we find that, due to the structure (41) of Vy, it can 
be S— deformed into qy given in (34)- (35), thus, for any normalized smooth test function of 
compact support in / [12] 



;$,h$) = 



\D$\ 



-dr + (k 



[n 



1)«) 



\<$>\ 2 H 
T r 2 



dr, 



(42) 



where = (-£-^ + S)<&, S = —fd\n(Ky)/dr. E in (40) is greater than or equal to a lower 
bound of (42). Since neither H nor D depend on ky, it is clear that (42) can be made 
negative for sufficiently high ky unless H is positive definite on I. In fact, if H < on 
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an open set O G I, the second integral in (42) will be negative for a test function with 
support contained in O, and ($,7Y$) < for high enough ky. On the other hand, if H is 
nonnegative in J, then E > since both integrals in (42) will be positive and 

*£-(n-l)K>0, (43) 
as can be seen by integrating by parts 

< / (VV J + V j V l )(V i V j + VjVi)y/§<Px (44) 

on the Riemannian compact manifold We conclude that a space-time is stable under 
vector perturbations if and only if H in (35) is non negative on I. As a first application, 
note from (35) that H — 1 in Einstein theory (a — 0), an already known [10] result on the 
stability of (1) under vector perturbations in General Relativity. Let us now consider the 
theories for which P(ip) in (5) has two real roots A 1 < A 2 , so that r extends to infinity (5). 
As in [12], define 

^Ai+As, A ^A i _A 2> x=t -± (45) 

in terms of which 

n — 3 n + 1 

~ 2(n-l) + 2(n-l)x 2 ' { } 

Equation (46) shows that all solutions are stable for these theories (note that that the Gauss- 
Bonnet term G( 2 ) b a in (2) is non trivial starting n = 3, so that (46) is positive definite in the 
relevant cases). If P(ip) has complex roots then A in (45) is purely imaginary and x 2 < 
in (46). For these theories the space-times (1) have a naked singularity [12]. The stability 
of space-times with a naked singularity in EGB theories is currently being studied [19]. 



IV. SCALAR PERTURBATIONS IN EGB GRAVITY 



After interpolating the n dependence in the perturbative equations for n < 9, we obtained 
a set of equations equivalent to 5G a i> = 0, in terms of the functions introduced in (21) and 
(26). The SGij = equations imply the condition 

= (F tt - F rr ) ((n - 2) (n - 3) a (/ - k) - r 2 + (n - 2) a rf) 

-2 (ra-2) (-r 2 + (n - 3) (n - 4) a (/ - k) + 2 (n - 3) a rf + a r 2 f") H L . (47) 

SG rt = is equivalent to 



= nr (2f-rf)^-f 



pjl TT 

rr , i 2 ur rt 2 U L 



dF rr 2 9F lv 
nr^r — h k s — — 2 nr 



5Gft = is equivalent to 

= [r 2 + (n - 1) (n - 2) a (« - /)] 



9F, 



dt 



drdt 



(48) 



n/r (-^ + 2r^) - 2 („ - 1) H L (k s 2 - nn) 



+ nr (jjf) ( r f + 2 (n - 1) /) - F rr (n (n - 3) / + k s 2 + nrf) 

-F rr + 2 r-^- J (2 r - (n - 1) (n - 2) a f) 
+2 (n - 1) (n - 2) a H L (k s 2 - uk) (2 k - 2 / + rf ) . 



(49) 
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5G ir = 0, is equivalent to 
= 2rf 



8Hr 



2 (n - 1) (n - 2) a (2 k - 2/ + rf) + F„. (2 r - (ra - 1) (n - 2) af) 



[r 2 + (n - 1) (n - 2) a (k - /)] 



2r/ 



2 (n — 1) 



<9r 



<9r 



<9 2 F, 



+F rr (rf + 2 (n - 3) /) - 2 r-^f + F« (-2 / + r/') 



(50) 



= 0, is equivalent to 

= (( n - 1) (n - 2) a (k - /) + r 2 ) 



<9F rf . s „ „, s d 2 F rt 



-2 (n-l)r- 



9t 



dF„ 
~dt 



+ 2 (n - 1) (n - 2) r— —a (rf - 2 / + 2 k) 



+r/(2r -( n -l)(n-2)a/')^. 



(51) 



Finally, 8G rr = gives 

= (r 2 + (n - 1) (n - 2) a (k - /)) |2 



nr 



dF tt 



„ 2 d 2 H L 
dt 2 
d 2 F 



f 



nr V (rf + 2 (n " 3) /} 



~ nr -Q^f ~ nF rr ((n - 3) / + r/') + 2n^f r - 2 (n - 1) H L (k s 2 - n«) + A; s 2 F ft 

-F rr + 2 r-^ J (2 nr - n (n - 1) (ra - 2) a /') 

_ 2 (n - 1) (n - 2) a fH L (k s 2 - n«) (2 « - 2 / + r/') . (52) 

These equations are used to solve for the different variables, until we get a single equivalent 
differential equation on a function of r and t, suitable for the stability analysis and other 
applications discussed in the Conclusions. We first solve (47) for F tt (t,r), 



Fa — F r 



2 (n - 2) (-r 2 + a {{n - 3) (n - 4) (/ - «) + 2/V (n - 3) + r 2 /")) 

a (n-2)((n-3)(/-/t)+rf)-r 2 L 1 ' 

Integrating (48) with respect to t gives 

f nr or 

plus an arbitrary function of r that we absorb into F rt using the freedom of adding to it an 
arbitrary function of r (see (26)). Inserting 

F rt =: j($ + 2H L ) (55) 
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x 



in (54), and the resulting expression in (49), we get an equation that can be solved for 
in terms of $, and <9$/<9r. The result is 

H L = [(« (n - 1) (n - 2) (/ - «) - r 2 ) (n (-rf> + 2 /) - 2 A; 5 2 )] 1 

jn/r (a (n - 1) (n - 2) (/-«)- r 2 ) ^ 

+ [(„ _!)(„_ 2 ) (/n ((n - 3) (/ - «) + r/') + k s 2 (/ - «)) a 

-r 2 (fc s 2 + n/(ra- 1))] $} (56) 

At this point we have explicit expressions for F tt , F rr , F rt , and Hl, all in terms of $(t,r) 
and its first and second derivatives with respect to r. Replacing these expressions in (50) 
and (52) we end up with equations that contain <9 3 $/(<9t 2 <9r), and <9 3 $/<9r 3 . However, the 
linear combination 

(-l/2)wM + g 
nr 2 (n — 1) (n — 2) (/ (r) — k) a — nr 4 ' 

where A and £> are the RHSs of (50) and (52) respectively, eliminates both terms, giving an 
equation of the form, 

<9 2 $ „ <9 2 $ „ 9$ 



m> /V + Qs * +te * = a < 58) 

where Qs and q$ are functions of r. One can show that any solution of (58), when inserted 
back into the formulas for F tt , F rr , F rt and Hl in terms of $, solves all the 5G a i> = equations. 
Note that (58) is of the same form as (32) for the vector perturbations. However Qs and qs in 
(58) depend not only on r, but also on all the parameters of the theory in a rather complicated 
way. In particular, we do not find the "factorization" property for the dependence on ks 2 , 
neither we find that qs(r) has a definite sign. Therefore, the S-deformation method does not 
appear to be readily applicable for the stability analysis in this case. On the other hand, 
and perhaps a little surprisingly, we have found an explicit form for the transformation that 
puts (58) in a standard Regge-Wheeler-Zerilli form. Namely, we introduce a function 0(t, r) 
such that 

$(t,r) = K s (r)i(t,r). (59) 
Then, choosing the integration factor 

r -V2n+i (nrf -2fn + 2k s 2 ) , x 

K s r = V ; 60 

v / r 2_ {n _ 2)a [rf , _( n _3) («-/)] 

and switching to tortoise coordinate r* (37) cancels the terms in d(j)(t,r)/dr* and yields 
an equation of the Regge-Wheeler-Zerilli form, which, after separating variables <f>(t, r) = 
0(r)e w ' gives (as in (40)) a stationary Schrodinger equation 

H<P = + V s <j> = -oo 2 <P = E<p (61) 

Unfortunately, the explicit expression for the scalar "potential" Vg in terms of r* and the 
parameters of the static solution is extremely long and complicated, and we have not been 
able to put it in a form that would be useful for a general analysis of the stability problem. 
It should be clear, nevertheless, that, for any choice of parameters, including n, Vs can be 
straightforwardly recovered, e.g., by means of a symbolic manipulation program, following 
the procedure outlined above. In the next subsection, some examples are analyzed by 
assigning particular values to the parameters. 
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A. Application: a scalar instability in small mass, 5D spherical EGB black holes 

For A = and k — 1, the e = — 1 branch of (6), 



reduces to the n + 2 dimensional Schwarzschild-Tangherlini [14] (Einstein) black hole in the 
a — > limit: 

/(r) = 1 - + 0(a) (63) 

In this section we will apply our results to exhibit a low mass scalar instability of 5D 
Schwarzschild-Tangherlini-EGB black holes. 

Before that, we note that, in general, a (assumed strictly positive from now on) can be 
used to introduce dimensionless quantities 

a ■= V A := Aa x : = roT 1/2 (64) 

and that / in (6) depends on a only through x, \x and A. If we define 

r* f x dx' 
x * .— _ = / ux (55) 

J Xo f(x',v,A) 



we find that (61) is equivalent to 



dx* 



+ aV S (p = aE<f>. (66) 



Furthermore, it can be shown that aVs = V(x,fi,A,ks 2 ), so that the stability problem re- 
duces to determining if the a— independent potential V(x,fi, A, ks 2 ) admits negative energy 
eigenvalues. For n = 3 ST-EGB black holes 



x 2 



f(x,fi) = l + — 



3x 4 



(67) 



and we must assume jl > 3/2, since there is no horizon below this mass value. This is a 
special feature of n = 3 EGB, for higher dimensions there is always a horizon [12]. The 
integral (65) for / given in (67) can be solved in closed form in terms of hypergeometric 
functions, and values of x* and V can be obtained for different x's and used to generate 
parametric plots of V vs. x*. All graphs in figures 1-6 were generated setting x Q = 2xh in 
(65), so that x* = when x = 2xh, and x* — > — oo as x — > x H + ■ Since V(x*) is bounded 
and lim x »_ += |- 00 V(x*) = 0, a sufficient condition for the existence of a bound state of negative 
energy is [15] 

poo 

V dx* < 0. (68) 



The above condition is certainly met in the graphs shown of figures 1 and 2, which exhibit 
the £ = 2 and £ = 10 potentials for p, = 1.7 [20]. For higher mass values (68) may not be 
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satisfied, and still the potentials may allow negative energy bound states. This is illustrated 
in figures 3-5. For these examples, we have fit a quadratic curve around the potential 
minimum, and checked that the ground state gaussian wave function (j> of the associated 
harmonic oscillator (plotted together with aVs) has a negative energy expectation value 



This shows that the spectrum contains negative energy eigenfunctions and that the 
black hole is unstable. For higher mass values in n = 3 (i.e, spacetime dimension five), 
as well as for higher spacetime dimensions, we were not able to find test functions 
with negative "energy" expectation values, there seems to be no scalar instability in 
these cases (see, e.g., Figure 6), although a more systematic study has not been com- 
pleted yet. Note however that a low mass instability under tensor mode perturbations 
was previously found in six dimensional Schwarzschild-Tangherlini-EGB black holes [11, 12]. 



The study of the linear stability of static solutions to the Einstein-Gauss-Bonnet grav- 
ity with spatial slices of the form x M + , E™ an n— manifold of constant curvature k 
can be carried out using the techniques introduced by Kodama and collaborators. In the 
classification of Kodama, et. al. [10], general linear perturbations can be constructed from 
appropriate harmonic tensors on E™, and classified accordingly into tensor, vector and scalar 
modes. Tensor perturbations were analyzed in [11, 12], the other modes being the subject of 
this paper. We proved that, as happens in higher dimensional GR stability problems [10, 16], 
the perturbation equations can be reduced to a single stationary Schrodinger-like equation 
on a function of the radial coordinate, the potential being different for the scalar and vector 
cases. Finding the linearized EGB equation and the appropriate integrating factors leading 
to the equivalent Schrodinger problem is far more difficult in the vector and scalar modes 
than in the tensor modes. However, after performing the calculations in spacetime dimen- 
sions five to eleven, the dimension dependence of the equations was interpolated by formulas 
that we believe are valid beyond this range and that reproduce the expected a — > (GR) 
results. 

In the vector case, we were able to prove stability using an S-deformation argument. The 
same result is found in higher dimensional GR. In the scalar case, however, we found an 
instability that is absent in GR [10, 16]. Although the complexity of the potential prevents 
us from drawing general conclusions, the EGB analogous to 5D Schwarzschild black holes 
are shown to be unstable below a critical mass. A similar result was found in [11], where it 
was proved that low mass 6D spherical, asymptotically Euclidian black holes are unstable. 
These results are relevant to TeV scale quantum gravity scenarios, where those black holes 
are predicted to be produced in high energy collisions [17]. 

The potentials that we have obtained have a number of applications beyond the study of 
the stability of black holes and cosmological solutions. Among the most immediate ones are 
the analysis of black hole uniqueness, quasi normal modes [18], and the analysis of stability 
of naked singularities in EGB gravity. These topics are currently under study [19]. 




(69) 
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FIG. 2: The potential V = aVs as a function of x* = r* / 'y/a, for n = 3, A = 0, p, = \ija = 1.7. 
The scalar perturbation corresponds to the £ = 10 harmonic. 
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FIG. 3: The potential V = aVs and a (non normalized) gaussian test wave function as a function 
of x* = r* I \/a, for n = 3, A = 0, jl = fx/a = 3. The scalar perturbation corresponds to the £ = 2 
harmonic. The normalized test function gives < —d 2 /dx* 2 >= 0.12 and < aVs >~ —0.28. The 
expectation value of the "Hamiltonian" is negative for this test function, implying the existence of 
negative energy eigenvalues. 




-8-1 



FIG. 4: The potential V = aVs and a (non normalized) gaussian test wave function as a function 
of x* = r*/y/a, for n = 3, A = 0, jl = ji/a = 3. The scalar perturbation corresponds to the t = 10 
harmonic. The normalized test function gives < —d?/dx* 2 >= 0.73 and < aVs >— —6.11. The 
expectation value of the "Hamiltonian" is negative for this test function, implying the existence of 
negative energy eigenvalues. 
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x* 



FIG. 5: The potential V = aVs and a (non normalized) gaussian test wave function as a function 
of x* = r* /\fot, for n = 3, A = 0, fi = /x/a = 6. The scalar perturbation corresponds to the i = 10 
harmonic. The normalized test function gives < —d?/dx* 2 >= 0.15 and < aVs >— —0.24. The 
expectation value of the "Hamiltonian" is negative for this test function, implying the existence of 
negative energy eigenvalues. 




FIG. 6: The potential V = aVs as a function of x* = r*/y/a, for n = 3, A = 0, jl = ji/a = 10. 
The scalar perturbation corresponds to the 1 = 2 harmonic. aVs has a small negative tail -shown 
in the inset- which cannot accommodate a negative energy state. 



